FLRW cosmology in Weyl-Integrable Space-Time 
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We investigate the Weyl space-time extension of general relativity (GR) for studying the FLRW 
cosmology through focusing and defocusing of the geodesic congruences. We have derived the 
equations of evolution for expansion, shear and rotation in the Weyl space-time. In particular, we 
consider the Starobinsky modification, f{R) = R -\- PP? — 2A, of gravity in the Einstein-Palatini 
formalism, which turns out to reduce to the Weyl integrable space-time (WIST) with the Weyl 
vector being a gradient. The modified Raychaudhuri equation takes the form of the Hill- type 
equation which is then analysed to study the formation of the caustics. In this model, it is possible 
to have a Big Bang singularity free cyclic Universe but unfortunately the periodicity turns out to 
be extremely short. 



I. INTRODUCTION 

In order to solve the problems (viz. singularity, 
flatness, horizon) concerning the early universe in 
cosmology, the modification of the classical geometrical 
structure of General Relativity (GR) becomes pertinent 
and various such modifications are attempted time and 
again ■ Various inflationary cosmological models Q , 
scalar-tensor theories and the models resulting from the 
string theory Q such as dilaton gravity are considered 
to solve the above-mentioned problems which are still 
not satisfactorily answered. These modifications amount 
in general to introducing a scalar field for some models 
while for others they appeal to a unified description of 
the fundamental forces for understanding cosmology of 
the early universe Q. In the context of unification of 
the forces, the Weyl space-time through its non-metric 
connection with a vector field (W^) could facilitate 
unification of the two long range forces, electromagnetic 
and gravity. 

The field defines an enlarged structure of space- 
time namely the Weylian manifold. The Weyl extension 
is a natural generalized space-time structure if we employ 
light rays and freely falling observers instead of clocks 
Q, [13] as the basic tools for maping the space-time, the 
EPS axiomatization shown [101 that Weyl space-time 
is the natural structure if we use light rays and freely 
falling particles instead of clocks as basic tools. In this 
approach, the Riemannian geometry is replaced by a 
Weyl geometry as the structure of the space-time which 
enlarges the framework of the GR space-time and as 
a direct consequence, a purely geometric fundamental 
new field is introduced into the theory of gravitation. 
The characterisation of the space-time and geodesic 
congruences with Weyl geometry in the Einstein-Palatini 
formalism need to be considered carefully for the com- 
plete manifestation of the gravitational fields [ll|, [l2| . 
The formation of caustics O [13] in a given space-time is 
important to examine the singularities of the space-time 



foliations [3, [lH • The caustics is an undesirable feature 
as it indicates failure of proper evolution of physical 
fields in certain region of space-time [l^ [13 ■ A caustic 
is therefore a singularity of a geodesic congruence where 
the equation of motion becomes untenable at such 
points [l^. The study of the formation of caustics 
in the field configuration relative a given space-time 
background is therefore very important and the analysis 
of the Raychaudhuri equation [l^, [H, [l^ i'^ diverse 
various settings [13, [l^l plays a crucial role in probing 
the physical viability of a model in relation to the 
occurrence of the space-time singularities [3, [T7| . 

In the present paper, we study the formation of caus- 
tics by examining focusing and defocusing of geodesic 
congruences in the Weyl space-time which becomes 
WIST for f{R) gravity in the Einstein-Palatini formal- 
ism. It is then equivalent to a conformal transformation. 
One of the basic tools of analysis is however the Ray- 
chaudhuri equation. The paper is organized as follows: In 
the next section, we set up the Weyl space-time structure 
which is followed by the computation of the kinematic 
parameters namely expansion scalar, shear and rotation 
with (1 -1-3) space-time split in section III. In Sec IV, we 
set up the Einstein-Palatini formalism for f{R) gravity 
followed by the study of formation of caustics in various 
cases by using the Raychaudhuri equation which is re- 
casted as a Hill-type equation in section V. Finally, we 
end up with the conclusion. 



II. WEYL SPACE-TIME 

In the Weyl geometry, the connection is related to the 
metric by the Weyl vector (W^) and is defined by 



(1) 



where V is a covariant derivative in the Weyl manifold. 
This equation is invariant under the conformal transfor- 
mations (i.e. gauge transformation) (23j 
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Using equation ([T]), one can write the connection (F^^,) 
in the foUowing form 



The Lagrangian in the classical electromagnctism can 
therefore be written as a geometrical term Q^F^yF^^ = 
V^yV^'^). But as pointed out by Einstein, such descrip- 
tion will also lead to a second clock effect (or a twin 
paradox of the second kind) which may cause inconsis- 
tencies (iB-E^l unless the tensor F^^y vanishes identically; 
i.e. 



r^, = f^, + A^, (3) 

where r(^p is the Levi-Civita connection and A(,'p is an 
additional term which depends completely on the Weyl 
vector, the so-called deviation tensor. 

For the torsion free connection, it also follows from the 
equation ^ that 

A^. - - ^ {S^Wy + S^W^ - g^yWP) . (4) 

Because of this additional vector, the Riemann tensor 
constructed from this connection is not antisymmetric 
in the last two indices. Therefore the Einstein-Maxwell 
theory can be described by the second Ricci tensor R^^p'^, 
also known as the homothetic curvature tensor or the 
segmental curvature tensor (V) [IJ] as given below, 

- If^.- (5) 

Here Ffj^y is the electromagnetic field tensor for the 
Weyl vector which serves as the electromagnetic four 
potential. 



d^Wy - dyWf, = 0. (6) 

It can not describe the electromagnetic field anymore, 
but it has various interesting cosmological consequences 
which will be of our concern in the rest of the paper. 
The equation ([5]) has the obvious solution for as a 
gradient of a scalar, = —dpcj)) [2^ This space- 
time is known as Weyl Integrable Space-Time (WIST) or 
reducible Weylian space-time. 

We will consider this geometry for present work be- 
cause it is the only viable subclass of Weyl Space-Time 
[30l |. The structure is the same as that of the conformal 
transformation of the metric .g^,y, the metric which de- 
fines the Levi-Civita connexion. In fact, if we define the 
metric (/) in the form 

= e*g^y, (7) 

the connection F^^ is now the Levi-Civita connection 
of the conformal metric /^j,. This is why the Weyl- 
Integrable space-times are also referred as conformally- 
Riemannian [3l| . 

It is then easy to calculate the Riemann and the Ricci 
tensors corresponding to the conformal metric and are 
given as below, 



(8) 



Rp.u{I) = RtLu{g) - V^^0 - i^gpu^ti 



(9) 



where is the covariant derivative for the Levi- 
Civita connection F of the metric g^^. 



III. (1 + 3)-DECOMPOSITION AND THE 
KINEMATIC QUANTITIES 

The study of the evolution of geodesic congruences es- 
sentially hinges on the evolution of the deviation vector 
^'^ connecting the two neighbouring geodesies in the con- 
gruence. Let be a tangent vector to the congruence of 
timelike geodesies and so we write Cu^ = C^u = which 



means that w'^V^^^ = ^"V^u^'. This leads to 

= (10) 
and thus is orthogonal to m'' everywhere. 

Since the congruence is associated with the timelike 
vector field u'^ , one can decompose the space-time metric 
into longitudinal and transverse parts with the projection 
tensor h associated with the unit vector u in the following 
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form, 



(11) 



We now have the fully orthogonally projected covariant 
derivative D defined by 



whieh measures the failure of being parallely trans- 
ported along the congruence. The evolution tensor on the 
subspace orthogonal to u can be decomposed in its irre- 
ducible components (i.e. into trace, symmetric-tracefree 
and antisymmetric parts) as follows: 



(16) 



'V pi 



h 



(12) 

It is easy to show that D^hap = h^Wphap, which 
is a Weyl space-time defined by the projection vector of 
Wp in the 3-dimensional subspace. We can now project 
the tensor field Vu in the components orthogonal and 
parallel to u as follows: 



u^_,Uu + - ( UuWf^ -f UpUyW ) (13) 



where W = u^W^ ~ —u'^Vy 



—(j). The usual decom- 
position is recovered for Wp =0. It is easy to see that for 
a timelike geodesic {ifV^Ui, ~ where V is the covari- 
ant derivative relative to the Levi-Civita connection), we 
have 



(14) 



taking equation ([T3| to the form 



w 



(15) 



It is worth noting that the Weyl vector measures the 
acceleration [iiy) in the Weylian space-time. These terms 
would be absent only when the Weyl field is zero. 

In order to have the kinematics of the geodesic congru- 
ence, one introduces the evolution tensor Bpi, = D^Up 



where the trace 0, the symmetric part Upu and the anti- 
symmetric part uipn are the expansion scalar, shear tensor 
(or vorticity) and rotation tensor respectively. This is the 
same as for the Levi-Civita connection, 

^yUp = DyUp = ^Ohpi, + + LUpt,. (17) 

After some straightforward calculations, it follows easily 



D,,u„ = D,,u 



1 . 



-hp,W. 



(18) 



In view of the above decompositions, the three kinematic 
quantities in the two cases are thus related by 



= e--w 



(19) 

(20) 
(21) 



Note that it is only the expansion scalar which is sensi- 
tive to the Weyl scalar field while the other two shear and 
vorticity remain unaffected. However the evolution equa- 
tions would be affected because the equations for shear 
and rotation also involve the expansion. 
The equation (fT9)) reduces to the trivial solution in the 
case of a WIST. In fact, in a FLRW spacetime with the 
redefinition of the scale factor as — > e'^a'^, we have 
0^0 + 3<p/2. 

Now in order to derive the Raychaudhuri equation, we 
begin with the evolution equation for Bp^, itself. 



J 



1 fl / 1 \ W" 

u^V^Bp, = -Bp^g""'B„, + u^i^R,;" + UpU^V,^ ~ -h^hj^ V pW^, ~ -W^Wp - — {B^, 



Bp^u,) (22) 



where 



hJ'hfu^R^^^^u, = u^u.R^,^" + UpU^V,^. (23) 
It is worth noting that the equation ([^^ takes the usual 



GR form u'V^Bpi, = —Bp^g'^'^ B^i, -\- u^UaR^up 

absence of Weyl vector Wp. Since u'Bp^V^g^'^ = —9W 
and using the equation ([2^ . the evolution equations for 
the three kinematic parameters are given by 
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(24) 

(25) 
(26) 



where R^y[pcr) — Vfiugpa-- Note that the usual part i.e. as 
in GR and the Weyl affected part separate out neatly in 
the above equations. Recall that the evolution equation 
for expansion, and similarly for the shear and rotation, 
the Raychaudhuri equation for the Levi-Civita connec- 
tion reads as 

u'^W^e = -rj^ ^ ~ Rp^u'^u'' (27) 

For the dynamical evolution, the "true" variable for the 



expansion is 9 as defined in and not 6 and also the 
time variation is defined as S^^, = u'^V-ySfj^i, 7^ u^'S/ ^S^^^, 
where S* is a tensor field. Therefore the evolution of a 
model defined in a Riemannian manifold with the pres- 
ence of a scalar field is different from its evolution when 
the scalar field is included in the definition of the con- 
nection. 

In the case of the WIST when the Weyl vector is a 
gradient, the equations (|24l25l26p reduce to 



2 11 
2 



In the particular case of a Friedmann Universe (i.e. in ab- 
sence of shear and rotation), the Raychaudhuri equation 
(|24p for the expansion scalar assumes the simple form, 

^ + y - ^0 = (31) 

which indeed reduces to the second Friedmann equation 
with 9 = 9 ^ 3a/a (where a is scale factor) in the stan- 
dard cosmology in absence of the scalar field. We em- 
phasize that by considering the equations (|9ll9l27p . we 
recover the equation pip and the same would be true for 
the shear and rotation as well. 



IV. EINSTEIN-PALATINI FORMALISM 

In the preceding section, we have derived the equation 
for the geodesic congruences in the case of a Weyl 



(28) 

(29) 
(30) 

I 

Integrablc Space-Time. So far we have not specified 
the Weyl field (/> ^tnd in fact it could be left as a free 
function which defines a new structure of the space-time 
or alternatively it could be derived on variation of the 
action in the Einstein-Palatini formalism (EPF). In that 
case it is equivalent to fixing the gauge ([5]). 

In fact all the actions of the following form in the EPF 
S = J d*x^f{R,i,,X), (32) 
where 1/; is a scalar field and 

X = -^gf^-d^^d,^ (33) 

are equivalent to a gravitational theory in a WIST. In 
the EPF, we do not fix any metricity condition as given 
by equation ^ and therefore it is not necessary to add 
a constraint in the action as it is done in the variational 
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principle with constraints |32| . The equations of motion 
are given by 



It then readily determines the Weyl scalar as 



(34) 

(35) 
(36) 

(37) 



which seems to show that the variational principle has 
selected a WIST. However, this is not a WIST but a 
Ricmann space-time with an undetermined gauge in the 
vacuum. What it says is that when matter is introduced, 
the structure of the space-time can no longer remain Rie- 
mannian but it becomes Weylian. This is an interest- 
ing feature that emerges from the consideration of the 
Einstein-Palatini formalism. 



V. EVOLUTION OF THE RAYCHAUDHURI 
EQUATION IN /(i?)-GRAVITY 

Let us consider matter as a perfect fluid described by 
the energy-momentum tensor. 



Tf^y = (P + P)Uf,U^ + Pgf_iy. 



(38) 



We assume that the matter fields are minimally coupled 
to the gravitational field {g) and the conservation equa- 
tion with the Levi-Civita connection T is 



V'-T^u = 0, 

which leads to 

p^e{p + P) = 0. 
Using equation ([T9|) . it takes the form 



p + e{p + p) = -{p + p)i. 

For f{R) gravity, the equation (|4ip writes 



(39) 
(40) 

(41) 



2MF(1 + w) 



2MF - p{l + w){l-Zw) 

3(1 + w)w 2 



-I- 



2MF - p{l + w){l-3w) 



p^ =Q (42) 



= 1 ( f'(R) 



where w = P/p, F = f'{R) and M = ^ ^ f"(R) 
R is solution of the equation 



- R 



Also the Raychaudhuri equation reduces to 



2MF 



+ 



3 2MF - p{l + w){l~3w) 

^ _ P-.f/2 



2MF - p{l + w){l-3w) 



F 



(44) 



We will consider a particular model where f{R) = R + 
/3R^ — 2A which then readily determines the Weyl scalar 
as 

(j> = ln{a + 213 p), with a = 1 + 8^ A (45) 

where we assumed the presence of dust {w = Q). 
The conservation and Raychaudhuri equation (|42l44p 
read respectively as 



a + 2Pp 

dp ir- 

a — pp 

9'^a + 2l3p 
3 a — Bp 



(1 - 8pA)p - Pp^ - 2aA 
2(a + 2/3p) ■ 



(46) 
(47) 



The standard results follow for /3 = 0. 

We will consider only the physical case where (/3, A) are 
positive which implies the positivity of a. Therefore by 
redefining the variable 6 = 3j{p)^, where 7 is a positive 
function of p, one can rewrite the equation (|47p in the 
following Hill-type equation. 



G + 



{2~a)p~(3p^ 



2Aa 



6(a + 2/3p)7(/3p) 



-G = 



(48) 



The caustic formation would occur due to focusing 
when G = and G < at a finite time [l^l- It is 
clear from the above equation that it is the sign of 
{2 ~ a)p ~ 13 p^ — 2Aa which will play the critical role. 
The special case where the cosmological constant is null 
simplifies the problem. In that case, for /3p < 1 it would 
be the usual GR result of focusing to the formation of 
caustics ultimately leading to the singularity. On the 
other hand; if however j3p > 1, the caustics formation 
will depend on the initial conditions. The focusing is 
then only possible for initially 60 < 0^ < where 6q 
is some critical value, else there would be dcfocusing. 
While in GR it is always focussing in the absence of 
rotation. 

In the case where we have a positive cosmological 
constant, two domains come out. One is for /3A > 1/24 
for which we can have defocusing (depending on the 
initial conditions) for any density p. And the second 
domain is for /3A < 1/24 where defocusing takes place 
for I3p < i-8M-VT^24?X f^^. > l-8M+VT^24?X _ 

We recover the previous result when A = 0. 



f'{R)R-2f{R)^-pil-3w). (43) 



In the domains where we found focusing there is always 
a singularity in the geodesic congruence. However in the 
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regime where there is a possible dcfocusing, one can avoid 
the singularity. The issue of singularity is thus sensitive 
to the choice of initial conditions. 
The equations (UHl) and (|T7)) after integration leads to 



0^ = 3A 



(49) 



1 + 8/3A + 2fip 

where C is a constant of integration which coincides 
with the spatial curvature (k) of the FLRW, C ~ 
—^k/{py^a'^) and the density is defined as p — p^a^/a?'. 
Therefore depending on the initial conditions, or equiva- 
lently on the spatial curvature, the model can lead or not 
to defocusing in the domains previously derived. We em- 
phasize that there is also an other initial condition which 
has to be chosen as p[t = 0) or 9{t = 0). 
It is interesting to note that 9 = 9 ""^^^p is singular at 

j3p = a while 9 is not. 9 and 9 are equal for /3 = 0. It is 
important to remember that the function which governs 
the expansion is 9 and not 9 and the two are however 
related through the equation ([T9|) . In the EPF, both 
metric and connection, which will however include the 
Weyl field, are independent variables. 
The possible solutions of equation ([^^ under different 
conditions manifest themselves in terms of four different 
branches (see FiglT]) for the evolution of expansion scalar 
with density. In the special case where A = 0, we have 
9 = Q (rebound) only for C < which means that the 
bounce can only be for a closed Universe {k = 1). This 
solution corresponds to the branches (III) and (IV) as 
mentioned in the FiglT] Depending on the initial condi- 
tion, the system can have two different behaviors. Let 
us consider, the case of collapse where initially fHp < a 
and as it proceeds p increases and 9 decreases. After a 
finite time, the density of the system reaches the value 
p ~ a/f3 for which the curvature is finite. This point 
corresponds to an attractor. In the case of GR where 
P Q, j3p ~ a corresponds to p ~ oo. This Universe 
can not have a primordial bounce except if there is a dis- 
continuous transition to an other branch. Therefore the 
singularity in GR is replaced by a discontinuity in this 
model. For the other branch (IV) where we integrate 
from Pp > a, the system is continuous but as the Uni- 
verse collapses the density decreases, after the bounce we 
will have an expansion of the Universe with an increasing 
density. The solution corresponding to the branch (III) 
can not lead to defocusing. However the solution corre- 
sponding to the branch (IV) leads to defocusing but in 
a Universe where the density increases during expansion 
(see FigU) and it can not describe our Universe. In the 
case when A = only these two branches are present and 
the corresponding solutions arc not satisfactory form the 
view point of our Universe. Therefore for the consistency 
of the model, a cosmological constant need to be included 
in it which leads to two more solutions of equation 
presented as branch (I) and branch (II) in FigHJ 

The branch (I) is a very well known result that ex- 
ist even in the case where f3 = 0. It describes a closed 
Universe with a cosmological constant. 



(I) 



(11) 



^ 





Figure 1: Identification of different branches for a closed Uni- 
verse in p-9 parameter plane corresponding to different generic 
features of the solutions of equation (|49|) with and without a 
cosmological constant. 



However the branch (II) represents a cyclic Universe. 
Let us consider a collapse of a Universe where the initial 
density is considered in the branch (II) . As it proceeds the 
density increases (see Figl2]) and reaches the maximum. 
This would then lead to defocusing and the collapse will 
end into a rebound {9 = 0) rather than a caustics. This 
is an effect due to the modification together with the 
cosmological constant. A necessary condition for the ex- 
istence of the cyclic Universe is < /3A < 1/24, which 
implies the extreme importance of the two constants (i.e. 
/3 and A) in the existence of the cyclic Universe. After the 
rebound, the density decreases in an expanding Universe. 
At small densities, the model will have an other rebound 
which creates the cyclic Universe. This solution exist 
only for a closed Universe and has a maximum periodic- 
ity T cx 100a/5. However the variation of the density of 
matter for the same Universe would be Ap oc 0.5Mpi//3. 
By considering /3 oc Hq^, the Universe will have a pe- 
riodicity larger that the Hubble time but at the same 
tine, the density of matter will have an extremely small 
variation Ap/ po oc 1 where po is the density today. 



VI. CONCLUSION 

One of the natural generalizations of the Riemannian 
nature of space-time is the Weyl space-time. Ehlers et 
al. had already shown that the Weyl space-time 

is indeed the most general structure for the study of 
gravitational field as a geometric field. Such general- 
izations are quite pertinent in addressing the present 
day observational challenges posed in terms of the dark 
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100 ;8'/' 

t 

Figure 2: Evolution of 9 and p for the cyclic Universe corre- 
sponding to a periodic solution. 



modification in f{R) gravity in Einstcin-Palatini formal- 
ism (EPF). In EPF, the Weyl field turns out to be an 
undetermined gauge and the space-time turns toWeylian 
from the Riemannian as soon as the matter is added to it. 

The issue of caustics formation is very critical to the 
location of occurrence of singularity in certain region of 
space-time and we have performed a detailed analysis 
of this and have applied it, as a trivial example, to the 
f{R) gravity cosmological models. It is found that the 
model R + f3R^ - 2A clubbed with the Weyl field does in- 
deed avoid the big-bang singularity which is already true 
without the /3— correction. For the Starobinsky modified 
gravity along with the suitable initial conditions, it is pos- 
sible to have a singularity free cyclic Universe which has 
unfortunately a very short periodicity. It is a pity that 
the model is not cosmologically viable, however it offers 
an interesting example of a non-singular cyclic model. 
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